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INTRODUCTION 

This report describes work completed or in progress during the second year of the 
three year grant NAGW-3075, A Theoretical Analysis of Thermal Radiation Form Neutron 
Stars, to Principal Investigator James H. Applegate. The bulk of the grant funds go to the 
partial support of a postdoc, Dr. Dany Page. The work described below is a substantial 
portion of Dr. Page’s research activities fro the past year. This work was carried out under 
the supervision of Dr. Applegate. 


Work Completed 

The unambiguous detection of thermal radiation from the surface of a cooling neutron 
star was one of the most anxiously awaited results in neutron star physics. This particular 
Holy Grail was found by Halpern & Holt (1992), who used ROSAT to detect pulsed X-rays 
from the 7-ray source Geminga and demonstrate that it was a neutron star, probably a 
radio pulsar beamed away from us. Two other pulsars, PSR 0656+14 (Finley, Ogelman, 
& Kiziloglu 1992) and PSR 1055-52 (Ogelman & Finley 1993), have also been detected 
in thermal X-rays by ROSAT. These results caused us to turn our attention away from 
developing ever more sophisticated cooling codes to applying our models to understand 
the new data. 

At an age of ~ 3.4 x 10 5 years, Geminga is in the photon cooling era. We have shown 
that its surface temperature of 5.2 x 10 5 K can be explained within the contexts of both 
the slow (modified Urea only in the interior) and fast (direct Urea, pions, kaons, etc) 
cooling scenarios; Geminga is too old to distinguish between these possibilities. However, 
agreement between the predictions and observations in either scenario is only possible if 
baryons are paired (BCS pairing, a 5 in superconductors) within most, if not all, of the 
interior of the star. In the slow cooling scenario, the surface temperature is too high 
unless the specific heat of the interior is reduced by extensive baryon pairing. In the fast 
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cooling scenario, the surface temperature will be much too low unless the fast neutrino 
cooling is shut off by baryon pairing (see Page & Applegate 1992). In this case the pairing 
must extend throughout the entire interior, and involve every particle species that can 
participate in fast cooling. These results are the strongest case to date for the necessity 
of baryon pairing in the interior of a neutron star; neutron pairing in the inner crust is an 
essential ingredient in the vortex line depinning theory of pulsar glitches, but this theory 
makes no statements about the baryons in the interior of the star. 

We also comment on the surface temperature determinations for the pulsars PSR 
0656+14 and PSR 1055-52, both of which are also in the photon cooling era. If the fast 
cooling scenario is correct, these stars also require baryon pairing throughout their cores. 
However, the observational uncertainties are larger for these pulsars than for Geminga, 
and a very conservative interpretation of the data using slow cooling and no pairing is 
allowed. 

This work has been written up and is accepted for publication in The Astrophyscial 
Journal (Page 1994). Preliminary versions have already appeared (Page 1993a, b). 


WORK IN PROGRESS 

All of our neutron star cooling models to date have used the unmagnetized effec- 
tive temperature-interior temperature relation for the outer boundary condition. We are 
improving the models by using published magnetic envelope calculations and assumed ge- 
ometried for the surface magnetic field to determine local interior temperature-emitted 
flux relations for the surface of the star. The magnetic field makes the electron contribu- 
tion to the thermal conductivity (the most important factor in determining the interior 
temperature-emitted flux relation) highly anisotropic; heat flows much more easily along 
the field lines than across them. As a result, regions of the surface where the magnetic 
field is nearly normal to the surface will be hotter than average and regions where the field 



-4- 


is neaxly tangential will be colder than average. Our procedure will allow us to include the 
effect of these temperature variations in our light curves. 

The four neutron stars for which there is good evidence that ROSAT has detected 
thermal emission from the surface (Vela, Geminga, 1055-52, and 0656+14) all have pulse 
fractions between 10% and 20% in the 0.08-0.50 keV band, which is the energy band 
most sensitive to thermal emission. From the work done so far, we conclude that there 
must be large temperature differences on the surface to account for these pulse fractions. 
These pulse fractions may provide interesting information about the configuration of the 
surface magnetic field of the neutron star. For example, models with a simple dipole field 
cannot produce pulse fractions as large as observed if the variation in surface temperature 
is entirely due to anisotropic thermal conduction. 

We have used blackbody emission at the local surface temperature for the emergent 
spectrum in our calculations. We hope to drop this approximation in favor of spectra of 
magnetic atmospheres in the future, although this undertaking may be beyond the scope 
of the current three year project due to its complexity. 

References 
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SURFACE TEMPERATURE OF A MAGNETIZED NEUTRON STAR: 
LIGHT CURVES AND SPECTRA FOR GEMINGA 


Dany Page 

Department of Astronomy, Columbia University 
Instituto de Astronomia, Universidad Nacional Autonoma de Mexico 


L INTRODUCTION 

The high quality ROSAT data from the neaxby neutron stars Vela, PSR 0646+14, 
Geminga and PSR 1055-52 present a new chalenge for theorists to provide good 
models for their interpretation. The published preliminary analyses of these data have 
shown that a large part of the soft X-ray received by ROSAT most probably come 
from the neutron star surface (in the case of Vela the contribution of the surrounding 
synchroton nebula has to be first separated). 

The strong magnetic field present at the neutron star surface affects the charac- 
teristics of the emitted radiation, significantly changing the spectrum as discussed 
by Pavlov and Shibanov in this workshop. Some of the resulting magnetized spectra 
have been used in data analyses and have led to substantially different conclusions 
about the neutron star surface temperature when compared with analyses using sim- 
ple black body spectra. All these works have so far assumed a uniform temperature 
distribution on the neutron star surface, with the exception of the inclusion of the 
contribution of the hot polar caps. However, another important effect of the magnetic 
field is to introduce a strong anisotropy of the heat transport in the outermost layers 
below the photosphere, resulting in a non-uniform surface temperature distribution. 

We present here some preliminary results of our modelling of this non-uniform 
surface temperature distribution. Beside generating composite spectra (from the 
temperature distribution) and phase dependent spectra (from the star’s rotation) 
we also naturally have the capability to produce light curves, thus adding a new 
dimension to the interpretation of the data. The prospect of this approach is not 
only to obtain more reliable measurements of neutron star surface temperatures, but 
also to learn about the configuration of the surface magnetic field. Our results so far 
use local black body emission; we hope to use more realistic spectra in the future. 




2. SURFACE TEMPERATURE 


All neutron stars of interest have an isothermal interior surrounded by an en- 
velope where a strong temperature gradient is present. Isothermality is reached at 
densities lower than 10 lo < 7 m/cro 3 , corresponding to a depth of at most 100 meters. 
Gundmundsson et al. (1983) have shown that the relationship between the interior 
temerature T, and the uniform surface temperature T, in absence of a magnetic field 

Ts , 6 ~ Til ( 2 ( 1 ) 

(T n = T/lO n K). This relation is mostly determined by the electron conductivity at 
densities ranging from ~ 10 4 to 10 s gm/cm 3 depending on the surface temperature. 
In presence of a strong magnetic field H the electrons are still free to move parallelly 
to H but their motion in the perpendicular directions is strongly inhibited. Detailed 
calculations have shown that the electron thermal conductivity is slightly increased 
along the field and strongly suppressed perpendicularly to the field. Compared to 
the zero field result, for a given internal temperature, the surface temperature will be 
slightly higher in regions where the field is almost radial but much lower in regions 
where the field is tangent to the surface. 

Due to the thinness of the layer where the magnetic effects are important the heat 
transport is still a one dimensional problem, at least in a first approximation. For 
a given surface magnetic field configuration H(0, <f>) the surface temperature T 3 (B, (j>) 
depends only on the angle 0 between the normal n and H as well as the internal 
uniform temperature X 1 , 

T,(M)=T,(0, Ti) (2) 

For parallel transport, 0 = 0, we use the results of Hernquist (1985) and for orthog- 
onal transport, 0 = 90°, we use Schaaf (1990a). For an arbitrary 0 Schaaf (1990b) 
has presented a two dimensional calculation and written his results in the form 

r,(0) = f(0) x T,(0 = 0) (3) 

where 

£(©) =£o + (l- fo) cos(0)°. (4) 

Thus £(0 = 0) = 1 and £ 0 = T s (0 = 90 °)/T s (Q = 0) < 1. The exponent a ~ 1 is 
itself a weak function of H and 7 1 ;. 
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3. GENERATING SPECTRA AND LIGHT CURVES 


Th procedure we use is as follows: 

• Choose the star’s mass M and radius R. 

• Choose a surface magnetic field configuration H(0, 4>). So fax we have only 
considered dipolar fields but more general configurations can easily be incorpo- 
rated. 

• Choose the internal temperature T,\ With this we generate the surface temper- 
ature distribution T s (9 , <f>) as described in Section 2 (see Figure 1). 

• Choose the observer’s position, i.e. the angle a between the rotation axis and 
the observer direction. 

• Rotate the star. We take 60 positions and at each position take a snap-shot 
of the star: we integrate the emitted flux over the visible part of the surface, 
taking into account gravitational lensing. This produces 60 phase dependent 
spectra which are then red-shifted. 

• Choose the stair’s distance D and the column density Npj for interstellar ab- 
sorption to obtain the phase dependent spectra as received at the entrance of 
the ROSAT telescope. 

• The spectra are passed through the PSPC resolution matrix to obtain the time 
(= phase) dependent predicted count rates. 

At each point of the stax’s surface we use black body emission at the corre- 
sponding temperature. Our spectra axe thus composite black body spectra. We hope 
to consider more realistic magnetic spectra in the future. 
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4. PRELIMINARY RESULTS FOR GEMINGA 

We borrow the basic characteristics of our analysis from the study of Halpern & 
Ruderman (1993, referred to as HR bellow), who considered both X-ray data from 
ROSAT and 7-ray data from GRO. The results we are presenting here have to be 
considered as illustrative: the number of degrees of freedom in our model is so large 
that we have not yet been able to do a comprehensive study. Moreover our method 
still uses local black body spectra. 

Beside the surface temperature distribution we also include two polar caps of 
adjustable size and temperature, located at the points where the magnetic dipole 
axis crosses the star’s surface. Within the restricted parameter search we have done, 
the configuration which gave the best fit to the published Geminga data (HR) is the 
following (see figures): 

• M = IAMq and R — 12 km (giving R°° = 14.8 km). No search was done 
on these two parameters, we just use typical values. The red-shift factor is 
e * = 0.8 and the maximum lensing angle 6 max = 120.7° (instead of 90° in flat 
space-time). 

• Off-center surface magnetic dipole field, in the equatorial plane (from HR: casi 
orthogonality of the dipole with the rotation axis is required by the 7-ray data 
within the outer-gap model; off-centering is needed to obtain a single peak in 
the upper energy X-ray light-curve). 

The center of the dipole is 2.4 km away from the star’s center and at an angle of 
60° with the radius (from fitting of the lower energy X-ray light-curve: distance 
from the center is restricted by the amplitude of light-curve, orientation with 
respect to the radius is used to reproduce the slight skewness of the light-curve) 

• Internal temperature Ti = 4 x 10 1 K (from fitting of the lower energy spectrum), 
giving en effective temperature at infinity = 4.47 x 10 5 /f and an average 
temperature (at infinity) T~ e = 4.18 x 10 s K. The maximum and minimum 
surface temperatures are respectively T™ ax = 5.36 x 10 5 AT and T~ n = 1.20 x 
10 5 AT. 

• Two polar caps of diameter 0.44 degree and temperature T™ = 3.2 x 10 6 K 
(HR: from fitting of the upper energy component of the spectrum). 

• Distance D = 120 pc and column density Nh = 3.8 x 10 2 ° cm ~ 2 (from fitting 
of the lower energy spectrum). 

As can be seen from the figures 2 and 3 the fit of the lower energy light-curve is 
quite good as well as the fit to the spectrum. 
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5. COMMENTS 


From the figures and experiences from the other trials we have done we can make 

the following comments: 

1. The light curves are quite smooth (of course), and gravitational lensing can 
only make them smoother. The sharp dips that seem to be present particu- 
larly in the lower energy curve could be evidence for some magnetospheric 
absorption. 

2. HR argued that the relative flatness of the medium energy light curve could 
come from a compensation between the lower and upper energy curves which 
are about 105° out of phase. If the upper energy curve is interpreted as thermal 
emission from the two hot polar caps then the contribution of these caps to 
the medium energy curve is totally negligeable: the flatness of this light curve 
must be explained by another mechanim. Anisotropic emission due to the 
magnetic field effects could be such a mechanims. 

3. With all surface dipolar magnetic field configurations we have tried it was im- 
possible to obtain a shift between the maxima of the lower energy light curve 
(due to emission from the main surface) and the upper energy light curve (due 
to polar cap emission): this could be an argument against a dipolar surface 
field (in case anybody needs such an argument). 

4. The PSPC resolution matrix we use is a more recent version than the one used 
by HR. It seems that the spectral fits need lower temperatures than in HR and 
also larger column densities. The value of Nh = 3.8 x 10 20 cm -2 used in the 
figures is much too large for a distance of 120 pc. If this is confirmed, it would 
be a strong argument for the inadequacy of the black body emission. A H 
or He atmosphere spectrum with low interstellar absorption should give results 
similar to a black body spectrum with strong interstellar absorption (but with 
a different temperature and star’s distance). 

REFERENCES 

Gudmundsson E. H., Pethick C. J. & Epstein R. I., 1983, Ap.J. 272 p.286. 

Halpern J. P. &; Ruderman M., 1993, submitted to Ap.J. 

Hernquist L., 1985, MNRAS 213 p.313. 
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FIGURE 1. Temperature distribution on the star’s surface, area preserv- 
ing projection. The values indicated along the isocurves are local tem- 
peratures, non red-shifted. The surface magnetic field is dipolar, in the 
equatorial plane, 20% off center and it makes an angle of 60° with the ra- 
dius. One hemisphere is warmer than the other, producing the skewness 
of the light curve in Fie. 3. 


FIGURE 2. Observed and theoretical spectra, in the ROSAT PSPC de- 
tector. The data are from Halpern Sc Ruderman (1993). The surface 
temperature distribution is shown in Fig. 1 and the observer is in the 
equatorial plane, at a distance of 120 pc (Nh = 3.8 x 10 20 cm -2 ). The up- 
per energy component (E > 0.5 keV) comes from the two hot polar caps 
(Tcap = 4 x 10 6 K). The separate contributions from the main surface and 
the polar caps are shown as dashed lines. These results are similar to 
the results of Halpern Sc Ruderman (1993). The dash-dot line is a black 
body at the effective temperature T e ff = 5.6 x 10 5 K and the dash-triple 
dot curve is a black body at the average temperature T ave = 5.23 x 10 5 K 
(non red-shifted values). 


FIGURE 3. Observed and theoretical light curves, in the ROSAT PSPC 
detector. The data are from Halpern Sc Ruderman (1993). The contin- 
uous curves are the total emission, main surface -+■ polar caps, while the 
dotted curves show the contribution from the main surface only. Notice 
that in the middle energy range the contribution from the polar caps is 
negligeable. In the upper energy band the contributions from each polar 
cap are shown as dash-dot and dash-tripple dot curves. 
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on Geminga. Section 3 describes the genera] physics of neutron star cooling relevent 2.2 Geminga’s temperature . 

to our present purpose. Section 4 discusses the various fast neutrino cooling scenarios 
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the upper value being probably an overestimate. We refer to Michel (1991) and Lyne (1992) and Shibanov et at. (1992) have partially extended these results including the 

k Graham-Smith (1990) for a discussion of the reliability of the spin-down age as an magnetic field effects. The general trend of these results, for H or He atmospheres, 

indicator of the true age. 1S that there is an excess emission in the Wien tail of the spectrum, compared to a 

blackbody, and the excess falls within the EINSTEIN and ROSAT detector ranges. 


This excess is reduced if metals are present (due to absorption edges) or when the is L ^ = 4? rR 2 <rT* oc T? 7 , where the effective temperature T e is converted into an 

effects of the magnetic field are taken into account. As a consequence, the use of internal temperature T, with the Ti — T e relationship calculated by Gudmundsson, 
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fermion process is about ( ksT/Ef ) 2 ~ 10 -6 time weaker than a two fermion process, extensive magnetic envelope calculations of Hernquist (1985) and Van Riper (1988) 

which is what the direct Urea - modified Urea emissivities show, Participation of a for transport parallel to H show that for a given inner temperature the surface 
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cooling cases are associated with quite different EOS at supranuclear density and we 
An important ingredient is the above mentioned 7i — T t relationship, Gudmunds- will discuss them separately, 

son, Pethick k Epstein (19S2, 1983) were the first to present a detailed study of it For modeling the slow cooling we consider five different EOS s from modern calcu- 

but did not include the effect of the magnetic field H which enhances the heat trans- lations, relativistic (MPA: Mutber, Prakash k Ainsworth 1987), nan-relativistic (FP. 

port parallel to H and strongly suppresses it in the perpendicular direction. The Friedman k Pandharipande 1981; WFF(avl4): Wiringa, Fiks k Fabrocini 1988) and 


two parametric EOSs PAL32 & PAL33 (PAL: Prakash, Ainsworth & Lattimer 19S8) two facts not taken into account in the models, and the stiffness of this EOS i: 

whose properties are intermediate to the above three. Each one of these calculations inseparable from to the lattice structure. 
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The PS EOS has already been considered in previous works (Nomoto & Tsuruta dramatic effect on the cooling, because it suppresses both the neutrino emission and 

1986, 1987; VanRiper 1991) and we will simply quote their results below. However, the specific heat. The pairing is in the >S 0 partial wave at low density and then shifts 

m this extremely stiff EOS the neutrons form a three dimensional lattice, thus having t0 the 3 P 2 partial wave at higher density. The protons in the core and the neutrons 

a totally different specific heat and different neutrino emissivity than liquid neutrons, in the j Mer crust ^ expected to be paired in the , 5o partia] wave whj]e the neutron 



pairing shifts to the 3 P 2 partial wave in the core. Pairing of protons in the 3 P 2 partial condensation give values of T c for 3 P 2 neutron pairing in between the two extremes 

wave seems to have never been considered, probably on the ground that its critical shown in Figure lb as T72 and T<2(m“ — 1), For 3 P 2 neutron pairing, background 
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free mass value. (The Hoffberg et at 1970 calculation used m* = 1). Of course m m 
is also poorly known at high density. Further extensive calculations within the ALS 
model (T&katsuka Tamagaki 1982 and references therein) with and without pion 


4. FAST NEUTRINO COOLING SCENARIOS would then be controlled by hyperon pairing. 

The models of Page k Baron (1990), Umeda et al (1993b) and Page k Applegate 
Fast neutrino cooling encompasses a variety of different scenarios, kaon or pion (1992) of kaon, pion and nucleon direct Urea cooling respectively with superfluidity 
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paired for the star not to drop into invisibility. In particular, if both A and E~ are multiplied by 0.1, are both acceptable within the observational uncertainty in the 

present together they undergo a purely hyperonic direct Urea process and thus one of Geminga data, 

them must be paired, One can expect that T c is lower for hyperons than for nucleons It should be mentioned that internal heating by friction of the crust neutron 

if the former have weaker interactions than the latter, and the star's temperature superfluid can significantly alter the thermal evolution of a neutron star when its 



temperature, thus its specific heat, is low. The models of fast coling with pions and in a 1.4 A/® star at a temperature T — 10 9 K for our five EOSs. Pairing will suppress 

reheating of Umeda et a!. (1993a), without core nucleon pairing, can produce a star C v exponentially when T <C Tc and the corresponding specific heat will practically 




(1979), Of critical importance here is the total specific heat of the star, which de- increases because the small 25%) decrease of the specific heat in the photon 

pends on the EOS and chemical composition (proton fraction). Pairing of nucleons is cooling era is not sufficient to compensate the significant reduction of the earlier core 

essential here because of its suppression of the specific heat. Table 3 shows the con- neutrino emission (during which only the very slow nn bremstrahlung is unaffected) 

tribution to the normal (i.e, without pairing) specific heat of the various components which gives a high temperature at the beginning of the photon cooling era. When all 






invisible after a couple of hundred of thousands of years. Consequently, any neutron assumed distance of 8 kpc which w'as later reduced by a factor 3 (Green & Gull 

star older than that with detectable surface thermal emission has an unpaired bary- 1982): with this new distance the resulting temperature would be low enough to 
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we still consider both types of neutrino cooling, letting observation be the ultimate lower temperatures and require less reheating. If both neutrons and protons are paired 

JU0 £ e ' within the whole core the theoretical temperature with slow neutrino cooling is much 

Since analyses of the ROSAT observations of PSR 0G56+14 and PSR 1055-52 have lower than the 7 x 10 s K reported. When comparing with the reheating models of 




Finley conclude that this temperature is compatible with the slow neutrino coolinc /v-krtiirw*’! ^ ♦ a l /ir + +■ i* \ « * , . 

S cooling) or too low (fast neutrino cooling) compared to the observed temperature of 

models, but only for models without pairing in the core. This temperature is too high Ceminga 

compared to the slow cooling when neutron core pairing is included, but reheating With fast neutrino cooling nucleon pairing is needed to stop the early cooling 

may explain the discrepancy and moreover fits with non blaclc-body spectra will give which, without this, would produce a star of temperature much lower than what 



is observed. If the fast neutrino emission is from hyperonic processes it is possible than 5 x 10 4 years where the slow cooling scenario predicts temperatures higher 

that the suppression we observe is due to hyperon superfluidity. It is not possible to than 0.9 to 1.1 x 10 6 K, depending on the exact age, while at later times both 
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reheating is needed, but several possible mechanisms have been proposed and may 
be able to provide sufficient reheating. 

The above discussion show’s that in order to distinguish clearly between the fast 
and slow neutrino cooling scenarios we need observations of neutron stars younger 


ACKNOWLEDGEMENTS 

Process Name Prongs Emissivity Q v 

I am grateful! to T. Ainsworth, J. H. Applegate, J, P. Halpern, M, Prakash and M. — ___ (erg/sec/cm 3 ) 





O0 


<M 



EOS Maximum mass star 1.4 A/q star EOS Core components Crust components 

. M maT (M s ) PmazUjrr^) Zp(%) IVJcm) x p (%) ■ n p e n e 


iO lO (D (O 0 

o o © o o 


CM 04 CO eo CO 

o> © 3 © 2 


U< ~ < 

^ Ph P-< 


o^ooo 

T-t © © 

c4 <N -2 



*a 

"o 


U3 

P 

bb 

.5 

"O 

*0 

q 

0 


^ . O) M 

5 < co CO 

JSo- J -1 
fci s < <: 

pL( P-t Ph 

£ 


05 

3 

oT 

8 

< 


to 

CO 

CM 

ei 

CM 


Cl 

< 

8 


2 2 


1 

p 

*5 

a 

1 

8 

P 

M 

«T 

P 

PL 

w 

w 

w 

i 

te 

2 

CL 

< 

2 




sf 

05 

05 

ci 

d 

d 

05 


CO 

05 

2 

2 

rH 



o 


’0 

a 

3 

S 

-a 

u 


CO 

00 

CM 


cl 

< 

<*£ 

05 


T3 

a 

p 

6 

'O 

p 

O 


cl 

03 

w 


E 

0) 

a 

*7? 

X 


8 

no 


W 


£ 

2 


jq 

PL 


O 

r— 

05 


45 

no 

2 

. 

5 : 

2 

a 

o 

w 

TO 

9 

jq 

u 

2 




U 

CO 

_T 

2 

2 

w 

2 

IS 

4 ) 

2 

c 3 

2 

e 

2 

5 

< 

2 

PL, 

a. 


£ 

9 

x 

2 


2 

2 

o 

W 

no 

P 


s 

Si 

•53 

< 

w 

a 

o 

on 

2 

1 

xj 

s 

cL 

2 

•8 

2 

^3 

2 

2 

00 

05 

y—i 

d 

9 

5 

2 


& 

E 

1> 

ffi 


to 


bO ^ 

h CM 


K JH 


5 

-*-» 

.S 


jq 

PL 


9Q 

O 

o 


>, 

j: 

cl 

13 

p 

CM 

CO 

05 

r-4 

< 

2 


£ 

co 

2 

d 

$-T 

1 

2 
J4 

8 

jq 

C5 

« 

s 

2 

d 

< 


►o 



•8 


V 

bO 

no 

$ 

3 


■R 

JO 

H 

*tn 

>> 


§ 

2 

PL 



OT 

4 

8 

«3 

>, 

a 

Ih 

PL 

bO 

'V 

PL, 

j 

u 

4> 

’C 

- 

PL 



05 


8 


05 

i — 1 

r2 

'a 

8 


2 

p 

cd 


*8 

<3 

S~i 

o 

e 

o 

p 

P 

Si 

cc 

lb 

a 

s 

-4J 

'w 

1 

g 

Vi 

< 


jq 

3 


a 

d 

u 

d 

i> 

cl 


B S 

\3 CM 


co a 

.H ^ 

flf v 

CO ^ 
2 £ 
6 £ 
ci ? 
i | 

1 I 
«» 

>-> CO 

2 £ 

o 

2 w 

-• i 


<p 

8 

> 


P 

PL, 


P-r 

2 

*j 

5 

CO 

1 

6 

d 

< 


2 


PL 

>* 


a 


to 

to 


CL 


> 

d 

2 

U 

1 

2 


9 

43 

2 

2 

p" 

a> 

jq 

3 

Q 

2 

d 

13 

a. 

6 

nj 

o 

2 

d 

c 

SC 

g 

PQ 

> 

6 


« i 


CL 

< 


c3 

O 

w 

2 

u 

cl 

w 

2 


a 

p 

6 

E 

o 

o 

<u 

> 

*C 

Pu 


►—3 

H 

2* 

♦j 

o 

if 

w 

a 

< 


CM 

8 

PQ 


rO 

PL, 

8 

05 


Oh 

43 


hQ 

H 

2 

i 

a 

< 


t- 

CO 

Tf 

< 


jq 

PL 

75 

P 

2 


CM 

< 


M 

PL, 

13 

p 

2 

to 

CO 


•s " w 

5 0 0 

-19 9 


Q Q 

43 43 


a 

< 


p 

a> 

v> 

’O 

fl 

p 

a 

< 


CL 

6 

o 

-a 

H 

2 

d 


p 

is 

S 

d 

d 


p^ 

< 

2 

«T 

o 

-o 

'O 

a 

d 


P PU 

a, >, 


^ to 


p 

o 

2 

u 

•o 

p 

<! 


£ 

js 

p 


Cl, 

< 


I 


CL, 


J 5 

•a 

d 

a 

PQ 


»-5 

Cl 

< 


-o 

p 

<e 


p 

co 

«y 


a 

•8 

-*-» 

0) 

d 

d 

I 

PQ 


CL 

< 


w 

< 

3 

? 

J 3 


CO 

<3 


13 

PQ 


JH .> 


bO 

T3 


PQ 


^5 

P 

V 


< 

«y 

d 

<u 

to 

p 

’3 

PQ 


r H 
PQ 

2 ^ 


bO 

rs 

*c 

a 


o 

PQ 

d 

p 

■<y 


Cl 

W 

2 


CL 

2 


p 

"9 


o> 

PQ 


H 

PQ 


< 

2 

8 
9 
9 
o 

2 

> a 

<i •g” 

* 3 

bfl 

TS 


PL, 

d 

<lT 

PL, 

2 

9 

9 

H3 

o 

P 

PQ 

2 

d 

a 

5c 

o 


J3 

a, 


?, 

CM* 

f'- 

05 

rH 

2 

d 

bO 

a 

>< 


a 

< 


C*H 

2 

9 

rG 

o 

d 


bO 5, 
p 

t* ~r 

2 $5 

Oi 05 

N3 05 

r-H 

^ 2 

^ 2 


CM CO 

hQ < 


a 

< 

CM 

05 

05 

f— 4 

o 

p 

o 

*N 

3 

d 

2 


S -g 


rP 

CL -i 

^ X 


p 

P', 


np 

9 

PL, 


Ol 

8 


00 *2. 
05 » 

*— • <, 


2 

> 


> 

d 


PQ PQ 


O 

•-M 

a 

jq 

o 

bO 

o 



d 

CO 


2 

2 

d 

2 

2 

2 

*-> 

i 

2 


bo 

bb 

>n 

g 

9 ^ 

8 

p 

4 ) 

§ 

15 

*X5 

15 

B 

JO 

J 3 

jq 

H 

E 


O 

o 

u 

2 

2 

2 


Glendemung, N. K. 19S5, ApJ, 293, 470 Michel, F. C. 1991, Theory of Neutron Stax Magnetospheres, (Chicago: The Uni 

Gehrels, N. k Chen, W. 1993, Nature, 361, 706 sity of Chicago Press) 

Gnedin, 0. Y. k. Yakovlev, D. G. 1993, Pisma v Astron. Zh. (Sov. Astron. Lett.), Miller, M. C. 1992, MNRAS, 255, 129 

19, (March 1993 issue) Muther, H., Prakash, M. k Ainsworth, T. L. 1987, Phys. Lett., B199, 469 



Niskanen, J. A. & Sauls, J. A. 1981, Preprint Schaaf, M. E. 1990a, A&A, 227, 61 

Nomoto, K. k Tsuruta, S. 1986, ApJ, 305, L19 Schaaf, M. E. 1990b, A&A, 235, 499 



00 

t— 

S 

t- 



s 

r-H 

lO 

1—1 

8 

d 

<£ 

■*T 

00 

o 

H 

d 

CO 


to 

bb 

o 


►o 

to 

CO 

K 

> 

a 

X 

>> 

Ah 

0 

1 

< 

cd 

Ah 

rG 

Ah 

cL 

00 

Ah 

oo 

cd 

rH 

§ 

8 

<7> 

W 



rQ 

Ah 

X 

H 

H 


> 






Ah 

W) 

u> 

rtj 

a 

1 

j* 

B 

e 

Ah 

s 

Ah 

S 5 

J 

(S3 


of 

r- 

CO 

I'- 


< 


a> 

r-H 

CD 

rH 




H 

H 

H 

> 


rtf 

i 

4f 


H 

CQ 


< 

.S 

1 — 1 


r 

W 


CO 

oo 

CD 

jd 

CO 

oo 

o 

§ 

Ah 

a 



CD 

CD 

*“f 

CD 

od 



< 


Ah 

< 

a 

o 

!"■* 

CD 

rH 

oi 

6 

£ 


2 

rtf 

H-* 

‘5. 

H 

rtf 

CO 

05 

CD 

rH 

HJ 

p 

V 

Q 

•i- 

Sf 

n 

2 

pi 


o 

Q 


cd 

rtf 

H-» 

rH 

§ 

1 

28 

bO 

a 

11,151 

rf 

Eh 

o 

to 

H 

Q 


I | : 

H H S 


Q IS P 

a £ s> 

«3 

Ah > CL, 


Prakash, M. 1993, in Proceeding of International Conference on ‘Realistic Nuclear Observation, ed. J. Ventura k D. Pines (Dordrecht: Kluwer Academic Publisher) 

Structure’, to be published in Phys. Rep. Wiringa, R. B., Fiks, V. k Fabrocini, A. 1988, Phys. Rev., C38, 1010 

Prakash, M., Ainsworth, T. L. k Lattimer, J. M, 1988, Phys. Rev. Lett., 61, 2518 
Prakash, M., Prakash, M,, Lattimer, J. k Pethick, C. J. 1992, ApJ., 390, L77 
Romani, R. W. 1987, ApJ, 313, 718 


jures cap ions paired according to Ainsworth, Wambach h Pines 1989. The temperature in ordinate 

is the effective temperature ‘at infinity’ i.e. red-shifted. The cross shows the Geminga 
Figure 1. a) Proton l S 0 pairing critical temperatures. CCY : Chao, Clark & temperature and age. 
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